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Transverse momentum dependent (TMD) parton distribution functions (PDFs), TMDs for short, 
are defined as the Fourier transform of matrix elements of nonlocal combinations of quark and 
gluon fields. The nonlocality is bridged by gauge links, which for TMDs have characteristic paths 
(future or past pointing), giving rise to a process dependence that breaks universality. It is possible, 
however, to construct sets of universal TMDs of which in a given process particular combinations 
are needed with calculable, process-dependent, coefficients. This occurs for both T-odd and T-even 
TMDs, including also the unpolarized quark and gluon TMDs. This extends the by now well-known 
example of T-odd TMDs that appear with opposite sign in single-spin azimuthal asymmetries in 
semi-inclusive deep inelastic scattering or in the Drell-Yan process. In this paper we analyze the cases 
where TMDs enter multiplied by products of two transverse momenta, which includes besides the 
PT-broadening observable, also instances with rank two structures. To experimentally demonstrate 
the process dependence of the latter cases requires measurements of second harmonic azimuthal 
asymmetries, while the py-broadening will require measurements of processes beyond semi-inclusive 
deep inelastic scattering or the Drell-Yan process. Furthermore, we propose specific quantities that 
will allow for theoretical studies of the process dependence of TMDs using lattice QCD calculations. 

PACS numbers: 12.38.-t, 13.85.Ni, 13.85.Qk 


I. INTRODUCTION 

When considering transverse momentum dependence of parton distribution functions (PDFs) one must account 
for process dependence, which is related to the color flow in the hard process and which is reflected in a process 
dependence of the gauge links (GLs) or Wilson lines that appear in the definition of the quark and gluon transverse 
momentum dependent (TMD) correlators [M3- In a field theoretical framework in terms of quark and gluon field 
operators, these definitions involve nonlocal combinations of such fields and hence necessarily also gauge links. To ‘feel’ 
the transverse momentum dependence they necessarily involve (covariant) derivatives or gluon field operators 052- 
Probed in the right way, however, through specific azimuthal asymmetries, often necessarily in combination with 
transverse spin asymmetries, the effects of transverse momentum dependence are not suppressed by the high-energy 
scale in the process pGlMTlbl . This involves processes that are sensitive to an observable transverse momentum, which 
introduces another scale into the process, allowing the unsuppressed appearance of operators that would be suppressed 
with inverse powers of the hard scale in inclusive deep inelastic scattering. Such processes should be based on TMD 
factorization 01 03 in order to link the observable to TMD correlators of quarks and gluons. 

In the parametrization of TMD quark and gluon correlators one encounters the TMD distribution functions depend¬ 
ing on x and p 2 . In addition transverse momenta may appear explicitly, e.g. in the simplest case there are terms in 
the correlator with a linear proportionality to the (relative) transverse momentum of quarks with respect to parent or 
produced hadrons. Examples of new terms in this case include single-spin asymmetries that at the level of distribution 
functions are hard to obtain without transverse momentum effects. They can be traced to particular time-reversal-odd 
(T-odd) matrix elements involving gluonic pole matrix elements [18H23J. Starting with GL-dependent TMD correla¬ 
tors, such matrix elements appear with calculable process-dependent gluonic pole factors depending on the path of 
the GL (34T3 l22l Iril - Eiil . The path can be future or past-pointing dependent on the color flow in the hard process. 
This leads for instance to the prediction of the sign change (factors ±1) of comparable single-spin asymmetries in 
semi-inclusive deep inelastic scattering (SIDIS) or the Drell-Yan process (DY) [Ml ■ 
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In this paper, we focus on the operator structure of the matrix elements in situations where transverse momenta 
show up at the quadratic level, either as p ^ or as a rank 2 tensor combination. The first one is relevant in studies 
involving the p% dependence of the TMDs and its possible process dependence. In the case of tensor combinations 
it leads to TMD functions in the parametrization that show up in the description of azimuthal asymmetries such as 
cos(2y>) or sin(2y;), where p is an appropriate azimuthal angle that can be defined in processes where at least two 
hadrons are involved, such as SIDIS or DY. 

The upshot of the paper is that transverse momentum dependence and the GL dependence in the description of 
TMDs in terms of matrix elements of nonlocal combinations of quark and gluon fields leads to a process dependence 
already for the rank zero TMDs, which are just the extension of the ‘standard’ collinear PDFs with an extra argument. 
In that case it affects in particular the p T -width of the TMDs. Process dependence of unpolarized TMDs has already 
been considered at small x in (27j-j3(|. In this paper, we outline the underlying operator structure, i.e. the split into 
universal operator combinations and its process-dependent coefficients. Besides the p T -width of TMDs also situations 
are considered in which higher-rank effects in p T show up, e.g. those involving the pretzelocity functions for transversely 
polarized quarks in a transversely polarized nucleon or in the description of linearly polarized gluons in an unpolarized 
nucleon. We link appropriate operators to a universal set of transverse momentum dependent distribution functions. 

In the analysis of operator structures it is convenient to use p T moments of TMDs, which formally are ill-defined. The 
average partonic transverse momentum (p^) is a theoretical quantity that beyond tree level requires a regularization 
and prescription in order to be well-defined. It is not directly observable. In this article we will not focus on this 
aspect and just view such transverse moments as idealized quantities, that are limits of well-defined Bessel moments 
for example Such Bessel moments are less sensitive to the large p T behavior of the TMDs, but they are still 

scale dependent and under changes in scale they will involve operator mixing between quark and gluon operators, and 
even between operator combinations describing unpolarized and polarized gluons. We will come back to the latter at 
the end of section UY1 

Assuming a proper definition in this way, and assuming TMD factorization, the p T moments can be connected to 
experimental observables. Of course, in practice the isolation of the average partonic transverse momentum from 
experimental measurements is extremely challenging due to the multitude of contributing effects. For example, the 
transverse momentum imbalance of dijet, dimuon and diphoton pairs in hadron-hadron collisions has been used to 
extract the average partonic transverse momentum [331 . 134j . Dimuons and diphotons are primarily sensitive to quark 
contributions, but under scale changes also gluon contributions enter. Dijets probe the transverse momentum of both 
quarks and gluons already at the leading order. In the analysis the observed transverse momentum Q T of the pair 
was divided into three parts: (Qt)p air/2 = (?4)intrinsic + (Pr)soft + (Pt)nlo, based on the fact that the ‘intrinsic’ 
transverse momentum distribution of partons beyond tree level is broadened by hard and soft emissions. Therefore, 
the theoretical description will require inclusion of “Y-terms” to bridge the low and high parts in the observable 
transverse momentum Q T . To reduce the sensitivity to the large p T contributions, here one could also consider 
extracting the Bessel moments directly from experimental data. See Ref. (Hj for a first analysis of this kind. 

There are further complications in the experimental extractions. Although there is no contribution from fragmenta¬ 
tion for these particular observables, the experimental transverse momentum resolution of the pair is a clear limiting 
factor and generally is of the same order or larger than the average partonic transverse momentum. In principle there 
are also contributions from spin correlation effects that were not accounted for in the experimental extraction pffil |37| , 
and another effect not accounted for in the extractions, is the process dependent color flow factors arising from initial 
and/or final state interactions (ISI/FSI). In fact, due to the presence of both ISI and FSI in dijet production, TMD 
factorization of that process is not expected to be valid [38] . The extraction of the average partonic transverse mo¬ 
mentum from the dijet imbalance in hadron-hadron collisions is therefore questionable, unless one can demonstrate 
that the spin correlation effects and the ISI/FSI effects are negligibly small. 

Despite the many complications in going from experimental measurements to trustworthy extractions of the pt 
widths of TMDs, it seems useful to investigate the possible effects of the color flow, and to find ways to assess the 
importance of such effects. The aim of this article is to shed light on this aspect and to propose ways to investigate 
it further quantitatively using lattice calculations, as we will discuss in section fVl 
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II. DEFINITIONS AND PARAMETRIZATION OF TMD CORRELATORS 


The quark and gluon TMD correlators in terms of matrix elements of quark fields 0,i including the Wilson lines 


U needed for color gauge invariance of the TMD case 


are given by 


quai'K 

01 


$ [ ^ ] (x,p T ;n) = j e ip< ( p ,SW’j(Q) U [0 ^i(0\ p ,S) 0 „ =o , 

2xTn) = J d ^ n p T (P,S\ F""(0) U m F™® ^, 0] \P,S)\^ 


(1) 

( 2 ) 


(color summation or tracing implicit), where we use the Sudakov decomposition p^ = xP M + Pt + crn^ for the 
momentum of the produced quark or gluon. In the TMD case, there are for a spin 1/2 nucleon in general eight 
leading contributing terms in the parametrization of the TMD correlator [3 j IH1 • but as already explained in the 
Introduction we focus only on specific examples, namely unpolarized quarks in an unpolarized nucleon or transversely 
polarized quarks in a transversely polarized nucleon, 


$ 


[Ed 


t \ / A 1 d/ 2-; , ,[£/]/ 2 \ <A . 7-L [U]/ 2 \ Prap^T 7x^75 \ P 

(x,p T ;n) = \x,p T )+h[ \x,p T )^ b $ T + \x,p T ) -^2- \~2 


(3) 


The spin vector is parametrized as S M = SlP^ + S% + M 2 S L n M (with in our case S L not needed). The symmetric 
traceless second rank tensor is p^P = PtPt ~ \p%9t^■ I n this parametrization this traceless tensor is used. Often 
one also encounters the expression in which just the symmetric combination is used. In that case one has a 

trace term that in the above expression is absorbed into hi introducing the combination /i0 = /t-0 + The 

notation indicates weighting with powers of —p^./2M 2 = p^,/2M 2 , of which the general definition reads: 

( 2 \ n 

7 ^) f...(x,pl), (4) 


and which upon integration are referred to as transverse moments. For gluons one has for the leading correlator in 
a polarized nucleon in general eight GL-dependent functions. We limit ourselves to unpolarized nucleons with for 
gluons the parametrization [4f| 


2xT fll '(x,p T ) = -g^ v f? [U,U] (x,p 2 T/ 


ILV 


PT_ h ±g[U,U'] ( 2x 


(5) 


which for gluons is sufficient to illustrate the complications when a product of two transverse momenta is involved. 
The (formal) integrated correlators are given by 


*«(*) = / ^e^<P,5|^(0)^^(e)|P,5)|^ T=o , 

2xT^(x) = /^ (P,S\ F^( 0) [ft F™(£) [ft |P,ft. n=(Ur=0 


( 6 ) 

(7) 


where the GLs are reduced to a unique straight-line GL which runs from 0 to £ along n, thus removing the link 
dependence. The relevant parametrizations for the collinear quark and gluon correlators in the examples (Eqs.[3]and 
0 which we study in this paper are 


$(x) = ^h{x) + hi{x)"/ 5 $T^ij, (8) 

2xT^(x) = -g^f!(x). (9) 

Here we do not discuss any complicating issues associated with scale dependence of the distribution functions and 
with convergence of the p T integrals. These are addressed briefly in [V] and more extensively in [32j,|4j]]. In the next 

sections, we make the GL dependence in the TMD functions like f^\x,p 2 ) for quarks or ff U ' U \x,p 2 ) for gluons 
explicit. 
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III. QUARK CORRELATORS OF RANK TWO IN A NUCLEON TARGET 


In this section, we study the quark correlator in Eq.[3]and we will identify the p T dependence with specific operators. 
For this we notice that the p T -integrated collinear correlator, 

J d 2 p T <& [u] (x,p T ) = &(x), (10) 

is GL-independent. Likewise one can consider higher transverse moments. Using the starting operator expression in 
Eq. m one finds rank two collinear correlators [(.>• 25j i 


d 2 Pi- $ [U] {x,p T ) = $do(x) 


i[u] 


E C GG,c Y GG,c 




a/3 


( x )> 


(ii) 


C— 1 


which are symmetric but not traceless. Subtracting its trace we get 

a/3 


J d2px M2 ® [U] ( x ’Pt) = ~^9t P *d-d( x ) 


2 

E 

C— 1 


rj\u] I g, a P ( T \ 

L 'GG,c \ V GG,c\ X ) 


^ Qt 3 &G-G,c( x ) 


( 12 ) 


These are relations involving on the right-hand side correlators with specific operator matrix elements, among them 
two (c = 1, 2) gluonic pole matrix elements multiplied with gluonic pole factors C^qq . The explicit operator 

structure of these matrix elements before doing the p T -integration is given by 




d^-Pd 2 ^ 




£-n—0 


(13) 


where the 0(£) operators are rank two combinations of id T (£) = iD“(£) — A“(£) and G“(£), defined in a color gauge 
invariant way (thus including GLs), 


K (0 = 9 / dp-P e(£.P - p-P) t/W] F™ {p)U^\ ,, (14) 

i r°° 

G“(C) = -/ drpP (15) 

" J — OO 


with e(C) being the sign function. Note that G Q (£) = G Q (£ T ) does not depend on £-P, implying in momentum space 
p-n = p + = 0, hence the name gluonic pole matrix elements. To be precise the correlator &gg,c contains the operator 
G“G^/M 2 . Including the quark fields in the correlator, the different color possibilities correspond to the color singlet 
combinations Tr(ipif>GG) (c = 1) or Tr('ip'ip)Tr(GG) (c = 2). The other correlator in Eq.[Tl] &dd, contains the operator 
(*9“)(«3 t)/M 2 . The collinear correlators and &gg,c are universal. In our particular examples, both of them are 
rank 2 symmetric correlators, but in general these are not traceless. 

In the next step, we are going to study the parametrization of the correlators. For this we will distinguish among 
the terms in the parametrization with different numbers of transverse momenta, i.e. terms with definite rank. To 
illustrate the procedure, it is sufficient to distinguish in our example the symmetric terms up to rank two, 


* [U] (x,P t) = $ [U] (x,pl) + $ a WV(x,pl) + ^ i > am (x,p* T ), 


(16) 


parametrized as 


$ [c/] (a;,^) + $“ [[/I(1) (x,p?,) = 


f[ U] (x,p 2 T ) 


h™ 


{x,pl)j5$ T 


t 
2 ’ 



(17) 

(18) 


where the tensor 4> Q/3 is symmetric but not traceless. Note that the weighted trace term [U( 1 )( a .^ 2 ) con _ 
tributes to the parametrization of the rank 0 part, while the trace term also has to be subtracted in the rank 2 part. 
This suggests absorption of the first two terms in Eq. [T7] into one term, which indeed is 

possible if we realize that the terms are GL-dependent anyway. We note that the notation using tildes above the 
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quark and gluon correlators in this expansion sometimes slightly differs from that in previous studies. In this paper, 
we simply use tildes for particular definite rank parts in the parametrization, which thus only depend on x and p 
The p T -integrated and weighted collinear correlators in Eqs. [10] HD and [12] are independent of the GL. The GL 
dependence is contained in the gluonic pole factors. It is then natural to also split the rank 2 contributions in Eq. [Tn] 
and start with the expansion 


<& [u] (x,p T ) 


H x ,pI) + ^J^®dd( X ’Pl) 


/ , °GG,c 
c= 1 


$ 


( 1 ) 

GG,o 


(x,p T ) 


PtolP 

~hp ®GG,c 


?( x i p T ) 


(19) 


Also in this expansion we thus must account for the trace terms appearing as (universal) matrix elements ^>g.g(x,p^,) 
and &g-g(xtPt) ' n the rank 0 part. In the remaining part indicated by the dots only terms with additional G-G 
will appear since, for the p T dependence, rank-2 structures are the highest ones for quark matrix elements in a spin 
1/2 hadron. Because $(x) and ®g.g(x) are GL-independent we can make the substitution <1>(x,p|) + &g[g (x,p%,) —> 
$(a;,p^). The rank 0 trace part of the gluonic pole correlator, however, must be kept and parametrized in analogy 
to the rank 0 part >1>(a;,p^). It cannot contribute to the integrated correlators, thus requiring that upon integration 
^G G d x ) = Following the naming convention in Ref. 9f] for the universal functions getting a superscript ( A ) for 

functions in $>gg and a superscript (Be) for functions in $>gg,c, this leads to the following parametrization in terms 
of universal functions, 



( 20 ) 

( 21 ) 


( 22 ) 


with 8f[ Bc \x ) = Sh[ Bc \x ) = 0. Making the trace term $g.g explicit by introducing additional functions Sf[ A ^ and 
Sh^ would be overcomplete 1 . 

Returning to our starting GL-dependent parametrization in Eq. [3] we have now constructed the universal set of 
functions, including fi and hiT, their p T -width effects 5f[ Bc ^ and Sh[ Bc ^ and three pretzelocity functions that give 
the GL-dependent distribution functions, 


fr = /i+E4lr ) +- 


AU] 


C—1 
2 


= h 1 + Y. C G ] Gjh[ Bc) 


C— 1 


MU] 

1 1T 




[U] h l.(Be) 
IT 


(23) 

(24) 

(25) 


Note that all functions in these equation depend on x and p For the pretzelocity distributions, this has been 
discussed in Q, but the broadening effects parametrized with two functions Sf[ Bc ^ is new. In principle traces of 
higher-rank tensors will give broadening effects coming with color factors C^ GGGG c , etc. These will not lead to new 
functions, but just to additional process-dependent broadening effects SSfi, SShi and also broadening effects Shf^ 


Note the missing factor of 1/2 for the terms in the last line of Eq. 29 in Ref. [{j. Also the order of some Dirac matrices in the Eqs. 31-33 
of the same reference should be reversed. 


l 
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in the pretzelocity functions. These satisfy 55fi(x) = 55 = 0, 55hi(x) = 56h^ = 0 and 5h^\x) = 0, etc. 

Such functions, however, only start playing a role in processes with rather complex color flow. Furthermore they just 
represent modulations in the p%, dependence. Comparing 55f\ with 5fi and fy, respectively, such modulations are 
expected to correspond to short distance effects in impact parameter space in the matrix elements. We expect these 
to become smaller and hope that this can be studied in lattice studies as discussed in Section [V] 


IV. GLUON CORRELATORS OF RANK 2 IN A NUCLEON TARGET 

As the second example, we consider the gluon correlator in an unpolarized nucleon with the GL-dependent 
parametrization given in Eq. [5] As for quarks the integrated and weighted results for gluon correlators give two 
types of matrix elements, 


d 2 PT r^ u ’ u '\x,p T ) = r^(x), 

' d 2 p T ^ v^ u ’ u \x,p T ) = r^f d (x) + £(7® r'^Oc). 


(26) 

(27) 


C—1 


For the gluonic pole correlators the index c now runs over 4 possibilities involving different color tracings. Includ¬ 
ing the two gluon fields in the correlator T (denoted with F) the relevant color structures are Tr(F[G, [G, F]]), 
Tr(F{G, {G, F}}), Tr(FF)Tr(GG) and Tr(FG)Tr(FG), respectively. These transverse moments can be used to iden¬ 
tify universal functions after introducing a parametrization in terms of definite rank, 


r^ u ’ u ']( x ,p T ) 


+ ^ r “Va‘W) 

4 

I V \v^ V t 1 ) (rr <n 2 \ _l_ P T( *P Ot£ 

I / . ^ dd n ^ d-d r\'* J ')P r r) I n rO ^ dC 


(28) 


with r (x,p%) absorbed in T^fya;,p^,), just as we did for the quark correlator in Ea. 1191 In the parametrization of 
the gluonic pole term a trace term T 9V G ) G c (x,p ^) is needed which must satisfy F^ G ) G c (x) = 0. With this expansion 
we can express the correlators in a universal set of TMDs depending on x and pfy 


r^(x,p T ) = -gr f?(x,p T ), 

r'GkV.Pr) = -«T 


/3}v aB uis 

Vt _9t_ 9t h ±g(A), 2 \ 


T^ c {x,p 2 T ) + tgfT 


a/3 -pnv 


2 VT - G G,c 


u{a B\v a/3 uv 

2\ _ 9t 9t ~ 9t 9t ,_ L g(Bc) / J2' 


( X,Pt) = 




(29) 

(30) 

(31) 

(32) 


where the pfyintegrated function 5ff(x) = 0. Returning to the GL-dependent functions in Eq. [5] we find 


/r ,, =/?+E4^r ) +- - 

c= 1 

ht 9[u ’ u> 1 = hpw+hp (Bc) +.... 


C— 1 


(33) 

(34) 


We already commented in the Introduction on the fact that the UV behavior of the functions requires care including 
operator mixture. The function h/ 9 is also a good example of the situation in which the large p T behavior of the 
function involves convolutions with functions corresponding to a different operator structure, that is in this context 
operators of a different rank. In this case it is the unpolarized gluon distribution that in the evolution equation also 
contributes to the large p T behavior of h/ a as already pointed out in Ref. |42, 43]. 
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V. GAUGE LINK DEPENDENCE IN LATTICE STUDIES 

TMDs are GL-dependent, but as the gauge links can be calculated for a particular hard process, one can express 
the GL dependence in gluonic pole factors. This will affect all TMDs, even the unpolarized quark and gluon TMDs 
causing a process dependence in their p 2 behavior and leading to nonuniversal T-even functions. To probe this process 
dependence in case of the unpolarized quark TMDs, one has to go beyond the comparison between relatively simple 
processes like SIDIS and DY, since the gluonic pole factors are both unity for simple future- and past-pointing gauge 
links, 


f[ + \x,p 2 T ) = f[ ] (x,p 2 T ) = f 1 (x,pl) + 8f[ m \x,pl). (35) 

In hadron-hadron scattering to hadronic final states one will find effects of unpolarized quarks coming from (among 
others) the correlators <f>l n +l and These gauge links are more complex than the simple staples t/J-^L which 

are consecutive Wilson lines 1/^ = 1/^1 ±00 -]Ujo T along minus direction or in the transverse direction 

at plus or minus light-like infinity. The GL t/P + l indicates a future-pointing GL that loops around once more, 
= U [^"^] U [^~q] U [^"^], while f7K D M indicates a future-pointing GL I/M and an additional traced GL t/^ D ^ = 

Tr({7[^]t/[c q])//V c [HUSH. We find in these cases 2 

f [ P + \x,p 2 T ) = fi{x,p 2 T ) + Q8f[ B1 \x,pl), (36) 

/j (D)+I (ir,p 2 ) = fi(x,pl) +8f[ m \x,pl) +8f[ m \x,pl). (37) 

We emphasize once more that the collinear integrated functions are the same in all cases. It should be mentioned 
though that strictly speaking one should not consider the collinear parton distribution functions as integrals over 
TMDs, as the latter require regularization of rapidity or light-cone divergences, and thereby do not necessarily lead to 
unique answers for integrals [4l| . In general, higher transverse moments, including the average transverse momentum 
(p 2 ) = 2M 2 f[ 1 \x), will diverge too, simply because the perturbative power law tail of the TMDs will not fall off 
sufficiently fast. For the purpose of regularizing this kind of divergence, a generalization of the weighting with powers 
of transverse momentum was suggested in |3l}, the so-called Bessel weighting. Bessel moments can be given as 
derivatives of Fourier transformed TMDs f(x, b ^) in impact parameter space, 

f~W( x ,b 2 T )=n\(-^d b ^ f(x,b 2 T ). (38) 


In the limit bx —> 0 the conventional transverse moments are retrieved, including their divergences, 
weighting regularized version of the average transverse momentum is given by 32]: 


The Bessel 


( P 2 T ){x,b 2 T ) = 2M 2 f[ l) (x,b 2 T ) =4./ d\p T \^^J 1 {\bT\\p T \) /iOuPt)- (39) 

Upon taking Mcllin moments, this quantity can be evaluated on the lattice [45|, [if|, which suggests a lattice study 
of the GL dependence of (cc, b^). For this purpose one can for example consider regularized versions of ratios 

such as 


f[ 1)[a+ \x) 

f[ 1)l+ \x) 


Ri 


1 + iV 


/r (1)[(D)+1 (*) . , ^ 

/i (1)[+1 (*) l + ^i’ 


(40) 


where R c = Sf[ 1 ^ Bc \x)/f[ 1 \x). The Bessel-weighted generalizations of these ratios can be evaluated on the lattice, 
e.g. 


/1 11(1)[D+1 (b 2 T ; p, 0 (PM 0 , 0 t ) 7 + ^ b T)\P) 

/j 1](1)[+! (b 2 ;/x,C) (P|?(0,0 t ) 7+ t/fo+i] ^(0,6r)|P> 


2 Compared to Ref. Q there is a redefinition of gluonic pole matrix elements and gluonic pole coefficients, leaving the product of the two 
unchanged. 











where the first superscript [1] refers to the lowest Mellin moment n = 1 and (1) to the Bessel-moment given above. 
These ratios offer a way to quantify the importance of gauge loops and of the process dependence of fi. Although 
it is not expected that the ratios are scale-independent, some of the scale dependence (of both /i and Q may cancel 
in the ratio a But even information about them at some fixed scales would already be very interesting. Any 
deviation from unity indicates first of all the relevance of gluonic pole matrix elements and second, the effect of the 
flux through a Wilson loop. How the latter changes with b is especially interesting. If the effect of an additional 
winding in the Wilson line enclosing the entire flux through the plane does not affect the average transverse momentum 
squared significantly, then one can conclude that the process dependence is not important to take into account for 
the unpolarized T-even distributions. Extensions to other T-even and T-odd functions are straightforward. 


VI. DISCUSSION AND CONCLUSIONS 

Even if quark TMDs in hadron-hadron scattering processes involve a combination of correlators $P + J, and 

some others with more complex gauge link structures, the consequences for process dependence of the unpolarized 
quark TMD f[ U ^ will most likely be hard to investigate experimentally. This will be looked at more carefully in a 
future study, which also includes other TMDs such as the T-even and T-odd rank 1 functions. Lattice calculations, 
however, could offer possibilities to investigate the gauge link dependence. For the gluon case, the situation may 
experimentally be somewhat less complicated, while the lattice calculation will be much more demanding in that 
case. The unpolarized gluon TMD ff U ' U ^ depends on two links and again the two simplest processes are sensitive 
to the same function: production of a colorless final state in hadron-hadron scattering, such as in the case of the 
Drell-Yan process or Higgs production, probes while cc production in SIDIS probes [37)] ■ Like 

f [ + ^ ^, it turns out that = f 9 ^ + ’ + \ which follows from a P and T transformation. For T-odd TMDs this 

equality does not hold |44j], e.g. for the gluon Sivers function /^ ^ ^ f^ + ’ + \ Processes in which other f^ U,U ^ 

functions contribute are for instance pp —>• H jet X and pp —>■ 7 jet X. In pp ^ H jet X [49] the subprocesses qq —»• Hg, 
qg —>• Hq , and gg Hg contribute, but when the momentum fractions of both initial partons is sufficiently small, the 
latter dominates, offering the possibility of accessing a single ff U ' U , which can be evaluated following the procedures 
as outlined in Ref. Q. A simpler way to probe a different f^ u,u 1 is in pp — * yjetX, if one can select a kinematic 
region where the partonic subprocess qg —»■ 7 q dominates over qq —»• 717 . In this case one can access /f^ + ’ ^. Such a 
study has already been proposed in pA scattering in [13 HU . Large A and small x (high energy) help to select the 
gluon induced subprocesses and lead to simplifications regarding the study of the process dependence of the gluon 
TMD ff U ' U ^ [13, HI, Ho]• In several of the mentioned processes also hf 9 ^ U,U ^ can be accessed [49)453 with the same 

link structure as f^ U ’ U ^ ■ A study of hf 9 ^ U ’ U ^ in eA and pA collisions has been performed in [52]. 

The average transverse momentum (pj.) and its broadening A p^ = (j>t)a — ( P%)p with atomic number A has been 
studied extensively in the literature in terms of the collinear factorization approach at twist-4 [53l - l58| . The relevant 
twist-4 parton distribution functions [59U60| are (at tree level) related to the first transverse moment f^\x) for quarks 

and gluons, but also determine the large-p T perturbative tail of the TMDs /and f^ u,u \ The collinear twist-4 
functions involve only light-cone operators, all having the GLs with finite paths along the light-cone and are thus 
process-independent. In the twist-4 calculations the process dependence is attributed to the hard partonic scattering 
factors. This implies process-dependent relations between the TMDs f\ U ^ and f 9 ^ U ' U ^ and the twist-4 functions, in 
analogy to the relation between the Sivers and Qiu-Sterman functions jfj). See Ref. [57] for a detailed discussion of 
the process dependence of nuclear broadening in the twist-4 approach. 

To conclude, in this paper we have elucidated the operator structure of quark and gluon correlators relevant to situ¬ 
ations where transverse momenta show up at the quadratic level, including as examples the quadratic p T dependence 
but including in this both the pj, dependence of the functions as well as the rank 2 tensor combination multiply¬ 
ing these functions. From this analysis it becomes clear which parts of the p T dependence of TMDs give process 
dependence. This is achieved by splitting the TMDs into universal operator combinations and process-dependent 
coefficients. This is of relevance (among other cases) for the p T -width of TMDs which enters in p T -broadening ob¬ 
servables, for observables involving pretzelocity functions for transversely polarized quarks in a transversely polarized 
nucleon, and for observables sensitive to linearly polarized gluons in an unpolarized nucleon. For the unpolarized 
gluon TMD f 9 ^ U ’ U ^ we have discussed ways to experimentally test the nonuniversality, in both electron-proton and 
proton-proton collisions. It would be very interesting if it could be established experimentally that even such a T-even 
TMD shows process dependence due to the gauge link structure. As it may be challenging to achieve this goal, we 
propose a way to study the gauge link dependence of f[ U ^ quantitatively using lattice QCD computations. 
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